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Abstract 

A brief review about the chemical composition of cosmic rays in 
the energy range 10 15 ^ E 10 20 eV is given. While there is convinc- 
ing evidence for an increasingly heavier composition above the knee, 
no clear picture has emerged at the highest energies, yet. We discuss 
implications about the origin of cosmic rays and emphasize system- 
atic differences related to data analysis techniques and to the limited 
understanding of the air shower development. 

1 Introduction 

When being asked to formulate the most crucial, yet unanswered ques- 
tions about CRs, one likely would end up with a list as such: 

1. Where and what are the sources of high energy CRs? 

2. How does the transport from the source(s) to the solar system 
affect their properties? 

3. What is the origin of structures observed in their energy distri- 
butions? 

4. Is there an upper end to the energy spectrum? 
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The enormous interest particularly in the highest energy CRs originates 
also from the fact that they provide a direct link to the problems of 
understanding the evolution of the early universe and to physics beyond 
the standard model of particle physics in an energy range not accessible 
to man-made accelerators. 

The chemical composition of CRs measured as a function of their 
energy provides a key to the answers of above questions and large ef- 
forts have been undertaken to provide reliable experimental data for 
comparison with models of CR origin. In fact, direct CR measurements 
performed at energies up to energies of several GeV have provided a 
great deal of information about the source composition, propagation 
effects, the 'age' of CRs, etc. However, for EAS experiments the sit- 
uation has been proven much harder than expected. This is mostly 
because of the large fluctuations of EAS parameters and because of 
uncertainties in corresponding EAS simulations. Nevertheless, signif- 
icant progress has been made in recent years and we have reached a 
point where high quality EAS data and models of CR origin and accel- 
eration can be confronted to each other. The information extracted 
from such measurements may be complemented by high energy (TeV) 
gamma and neutrino radiation. The latter messengers point directly 
to the source but are only secondary probes since they originate from 
interactions of CRs with the local environment at the source. Under- 
standing the common picture of photons, neutrinos and high quality 
CR energy and composition data may finally yield a conclusive picture 
to the open CR-questions of the century. 

2 Experimental Results 

The energy spectrum of CRs follows a simple power law behaviour 
(dJ/dE oc E~ 7 ) over many orders of magnitude. There are only two 
well established structures in the CR spectrum; the knee at E ~ 3 
PeV where the index changes from 7 ~ 2.75 to 3.1 and the ankle at 
E ~ 5 EcV where the spectrum turns up again. The origin of the struc- 
tures, however, remains unknown as is the existence or non-existence of 
the Greisen-Zatsepin-Kuzmin cut-off. Since structures in a spectrum 
generally contain important information, many experiments have fo- 
cussed to study these energy ranges. However, due to the low flux 
of high energy CRs, all of these structures are accessible only to EAS 
experiments. The primary energy and composition is then inferred 
either from particle measurements at ground or from Cherenkov and 
fluorescence light observations, or from combinations of both. Major 
particle arrays in operation or still delivering newly analysed data in- 
clude: AGASA, CASA-mia, EAS-TOP, GRAPES, Haverah Park, HEGRA, 
KASCADE, Maket-Ani, msu, SPASE, and Tibet-As7. Eas experi- 



ments exploiting optical observations include: blanca, dice, Fly's 
Eye, HEGRA, HiRes, and tunka. A brief overview about most exper- 
iments and their observational techniques can be found in Refs. 0, ^]. 

Previous results about the CR composition were reported e.g. in 
H Qj. The traditional and perhaps most sensitive technique to infer 
the CR composition from EAS data is based on measurements of the 
electron (N e ) and muon numbers (N^) at ground. For consistency tests 
or for improving the sensitivity, N e and may be complemented by 
additional EAS parameters, such as the electron lateral distribution 
(ldf), muon production height, hadronic observables, particle arrival 
times, etc. 

Cherenkov and fluorescence observations are restricted to clear 
moonless night which imposes a maximum duty cycle of approx. 10%. 
In non-imaging Cherenkov measurements the composition is inferred 
from the steepness of the lateral Cherenkov light intensity at distances 
up to 120 m from the shower core. Instead of presenting the Cherenkov 
ldf for comparison with EAS simulations, many experiments have pre- 
ferred to introduce another intermediate step by inferring the maxi- 
mum of the longitudinal shower development, X max , by means of EAS 
simulations so that only in a second step these inferred X max values 
are compared with EAS simulations. Obviously, this introduces addi- 
tional systematic uncertainties and may result in inconsistencies when 
for example different hadronic interaction models are applied to the 
(Ch-light intensity «-> energy )- and (Ch ldf <-> X maa; )-relations on 
the one hand and to the X max distributions of EAS simulations em- 
ploying different hadronic interaction models on the other. Up to now, 
only fluorescence measurements allow direct measurements of X max 
from observations of the longitudinal shower profile. 

A compilation of experimental results on X max is shown in Fig. 1 
for a wide range of primary energies The data are compared to 
CORSIKA simulations using different hadronic interaction models. As 
can be seen from that figure, differences between the models are on the 
order of 25 g/cm 2 in the knee region and increase to about 40 g/cm 2 
at higher energies. Systematic uncertainties in the data are estimated 
to at least be of similar size. Again, besides the shifts seen in the simu- 
lated A maa; -lines, also the experimental A maa; -values of the Cherenkov 
data may shift when extracted using different interaction models. With 
these uncertainties kept in mind, all data seem to suggest a mixed com- 
position over the whole range of energies with no clear picture emerging 
about a change in the composition. Only the CASA-blanca || and 
SPASE- vulcan |t| data may indicate a lighter composition towards the 
knee with a turnover to a heavier one above the knee. The Cherenkov 
imaging experiment dice B appears to contradict that picture. How- 
ever, electronic saturation effects were recently discussed to have biased 
their data ||. The weak indication from light to heavy when crossing 
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Figure 1: Compilation of X max from different experiments compared with CORSIKA 
simulations using different hadronic interaction models (adopted from Ref. ||). 



the knee is followed by a trend back towards a mixed-light composition 
in the energy range 10 17 - 10 18 eV in the HiRes-Mia data 0. The 
original Fly's Eye data suggest a heavy dominated composition with a 
much weaker change to a lighter composition only around E ~ 10 19 eV 
Jo) . In addition to experimental effects large uncertainties are added 
by the poorly understood hadronic interaction models at these high 
energies. 

Data from EAS arrays are mostly presented in terms of the mean 
logarithmic mass, (In A) = J2i r i^ nj ^i- Despite the drawbacks of this 
observable |lj] we shall follow this approach to ease the comparison of 
data from different experiments. Figure 2 shows a compilation of sev- 
eral ground arrays. Data from the balloon borne experiments Jacee 
[fl2"[ and Runjob |ll| and results from a simplified leaky-box [Q and 
diffusion model calculation assuming a knee at a constant rigidity of 4 
PV are included for comparison. All EAS arrays observe a change 
towards a heavier composition across the knee. However, absolute dif- 
ferences between experiments up to (In A) ~ 1 are found in the transi- 
tion region. As comparisons between QGSJET and neXus indicate |Tq] , 
another overall shift by (In ^4} ~ 0.5 should be added from uncertainties 
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Figure 2: Mean logarithmic mass from direct CR measurements and air shower 
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in the interaction models. The EAS data presented in Fig. 2 were all 
analysed with reference to CORSIKA/qgsjet. Reasonably good agree- 
ment is seen also between EAS data and direct measurements from 
the Runjob Collaboration favouring an almost constant composition 
from approx. 10 13 to 10 15 eV. Jacee on the other hand suggests an 
increasingly heavier composition starting already at approx. 10 14 eV. 
The cause of this difference is not fully understood. Helium fluxes 
in Runjob data are lower by a factor of 2 as compared to Jacee 
and, perhaps more importantly, Jacee attributes sub-Fe elements to 
the Fe-group thereby increasing the (In A) values. The KASCADE data 
shown in Fig. 2 result from an unfolding procedure of the N e - and 
iV* r -size distributions J|| and perhaps provide the best clue to the 
knee puzzle. The preliminary data convincingly show knee structures 
in energy spectra selected by primary mass groups. Furthermore, their 
knee position changes towards higher energies with increasing mass of 
the primaries. Detailed examination of the individual knee positions 
suggest a constant rigidity for all components |^3|, ^ as is expected 
in astrophysical acceleration and transport models. 

The only composition analysis above 10 17 eV based on a ground 
array and using state of the art EAS simulations has recently been pre- 
sented by the Haverah Park group . Their analysis is based on the 
shape of the observed lateral distribution function derived from the 
water-Cherenkov signals. Assuming a constant bi-modal composition 
in the energy range 0.2-1.0 EeV best agreement to the data is found 
for (34 ±2)% protons and the rest iron. Systematic uncertainties in- 
troduced by the choice of hadronic models were estimated to approx. 



14%. This error is included in the data point shown in Fig. 2. Within 
error bars, this finding is in agreement to the KASCADE results at lower 
energies. Furthermore, there may a turnover back towards a medium- 
light composition at higher energies. Such a change is indeed reported 
by fluorescence measurements of the HiRes prototype working jointly 
with the MIA muon array for energy determination [To[ and has been 
reported to some extend also by the Fly's Eye collaboration jll] (see 
above). However, quantitatively there are still significant differences in 
the composition at energies above 10 17 eV so that the question appears 
far from being resolved. 



3 Final Remarks and Outlook 

A variety of experimental techniques is applied to infer the energy 
and mass of CRs from EAS. In the knee region, two approaches have 
been proven successful: extended particle arrays measuring electrons, 
muons, and high energy hadrons on the one side and non-imaging 
Cherenkov measurements on the other. At the highest energies, the 
Cherenkov technique becomes impractical but fluorescence observa- 
tions take over. The wide scattering of reconstructed data points in 
any of the approaches, even when employing the same hadronic inter- 
action model, may hint to hidden experimental systematics and/or to 
incomplete data analysis techniques in at least some cases. Indeed, 
the data analysis is complicated by many subtle details: the energy 
spectrum is very steep but (i) shower fluctuations are in general very 
large, (ii) they depend strongly on the (to be reconstructed) mass of 
the primary particle and, if not corrected for, (iii) their effect mimics 
a primary dependent increase of fluxes which in addition depends on 
the steepness of the underlying elemental spectra. As a consequence 
of the latter, the mass dependent changes of the spectral indices at the 
individual knee positions causes a bias towards a lighter composition 
above the knee. Whether that effect is responsible for the apparent dip 
towards light masses seen in knee region of the Cherenkov data (see 
Fig. 1) cannot be said at present. Furthermore, on top of these com- 
plications, experimental reconstruction accuracies generally dependent 
on the primary energy and mass and have to be carefully corrected for 
as well. Clearly, sophisticated data analysis techniques are required, 
such as unfolding procedures which do not make any assumption about 
the shape of the individual primary energy spectra or about their rela- 
tive abundances. Very promising results have recently been presented 
by the KASCADE collaboration and the techniques may be ap- 

plied also to other existing experimental data. The present limitation 
of such approaches is, besides the available Monte Carlo statistics, 
mostly given by the reliability of the simulated tails in the used EAS 



observables. 

The situation at higher energies becomes more complicated mostly 
because of increasing uncertainties of interaction models. Akeno, Fly's- 
Eye [p5[ , and HiRes-Mia measure a composition rich in Fe at 10 17 eV 
which may become lighter at higher energies. On the other hand recent 
analyses of Haverah Park data indicate a mostly light composition 
already above 2 • 10 17 eV |§. 

With these limitations in mind, what can be learned from the 
present EAS data? There appears convincing evidence for a change 
from a 'standard' light (p+He dominated) composition to a heavy one 
(Si or heavier) in the energy range from the knee to about 10 17 eV. 
According to preliminary kascade data, this change in composition 
is caused by subsequent 'breaks' in the energy spectra of the different 
elements. The data also provide support for the knee being an effect of 
constant rigidity at E/Z ~ 4-10 15 eV, such as is expected in astrophys- 
ical models for confinement and/or acceleration in magnetic fields. The 
iron knee would then be expected at E » 10 17 eV. It will be interesting 
to confront such data with detailed model predictions. If the knee in 
the spectra is mostly due to galactic modulation, the question arises 
about to which energies the Fe-group could dominate the all-particle 
flux and to where the extragalactic component would take over. As 
recently argued by Berezinsky et ai, this transition could occur at en- 
ergies even lower than 10 18 eV so that the composition naturally would 
become light again. 

To conclude, enormous progress has been made in recent years in 
measuring the CR composition more accurately than in earlier gener- 
ations of EAS experiments. This is particularly true for the energy 
range of the knee but new analyses have also been reported at higher 
energies. Besides further improving the data analysis techniques, a 
wealth of new data at high energies will become available in the very 
near future from KASCADE-GRANDE, HiRes, the Pierre Auger Obser- 
vatory, and in the more distant future possibly from space experiments 
like EUSO and owl providing the required input for solving the open 
questions about CRs . 
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